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Post-Newtonian theory of motion of celestial bodies and propagation of light was instrumental 
in conducting the critical experimental tests of general relativity and in building the astronomical 
ephemerides of celestial bodies in the solar system with an unparalleled precision. The cornerstone 
of the theory is the postulate that the solar system is gravitationally isolated from the rest of the 
universe and the background spacetime is asymptotically flat. The present article extends this 
theoretical concept and formulates the principles of celestial dynamics of particles and light moving 
in gravitational field of a localized astronomical system embedded to the expanding Friedmann- 
Lemaitre-Robertson- Walker (FLRW) universe. We formulate the precise mathematical concept 
of the Newtonian limit of Einstein's field equations in the conformally-flat FLRW spacetime and 
analyze the geodesic motion of massive particles and light in this limit. We prove that by doing 
conformal spacetime transformations, one can reduce the equations of motion of particles and light to 
the classical form of the Newtonian theory. However, the time arguments in the equations of motion 
O ■ of particles and light differ from each other in terms being proportional to the Hubble constant !K. 

This leads to the important conclusion that the equations of light propagation used currently by 
»' I \ Space Navigation Centers for fitting range and Doppler-tracking observations of celestial bodies are 

missing some terms of the cosmological origin that are proportional to the Hubble constant 'K. We 
also analyze the effect of the cosmological expansion on motion of electrons in atoms. We prove that 
the Hubble expansion does not affect the atomic frequencies and, hence, does not affect the atomic 
time scale used in creation of astronomical ephemerides. We derive the cosmological correction to 
the light travel time equation and argue that their measurement opens an exciting opportunity to 
determine the local value of the Hubble constant "K in the solar system independently of cosmological 
observations. 
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I praise you, Father, Lord of heaven and earth, 
because you have hidden these things from the 
wise and learned, and revealed them to little 
children. 

Luke, 10: 21 



I. INTRODUCTION 

Post-Newtonian celestial mechanics is an important branch of modern gravitational physics [isl . [65j including 
gravitational wave astronomy [66j . It is also a practical tool of modern appUed astronomy used for precise navigation 
of spacecrafts in deep space and for calculation of high-precise astronomical ephemerides of celestial bodies in the 
solar system (s^ . [STj . The mathematical formalism of the post-Newtonian approximations is getting progressively 
complicated as one goes from the Newtonian limit to higher orders [6, 55]. For this reason the theory has been 
developed under a basic assumption that the background spacetime is asymptotically flat. Mathematically, it means 
that the full spacetime metric, (the Greek sub-indices a,/3,j,. . . takes values 0,1,2,3), is decomposed around the 
background Minkowskian metric, f = diag(— 1, 1, 1, 1), into a linear combination 

5q/3 — iat3 + hajj , (1) 

where the perturbation, hap, is the post- Newtonian series with respect to the powers of 1/c, where c is the ultimate 
speed in nature (the speed of light in vacuum). Post- Newtonian approximations is the method to determine h^fj by 
solving Einstein's field equations with the tensor of energy-momentum of matter of a localized astronomical system 
Tq^ as a source of the field, by iterations starting from hap = in the expression for Tq^. The solution of the field 
equations and the equations of motion of the astronomical bodies are derived in some inertial coordinates r" = {ct, r} 
where t is the coordinate time, and r = {r*} are spatial coordinates. The post-Newtonian theory in asymptotically-fiat 
spacetime has a well-defined Newtonian limit determined by: 

(1) equation for the Newtonian potential, U — c^/ioo/2, 

where p = TqOj is the density of matter producing the gravitational field, 

(2) equation of motion for a massive particle 

r" = VU , (3) 

where V = {d/dx^} is the operator of gradient, r = r{t) is a radius-vector from the origin of the coordinates to 
the current position of the particle, and a prime denotes a total derivative with respect to time 

(3) equation of motion for light (massless particle) 

r" = , (4) 

which tells us that light moves along a straight line with constant velocity, r — ^ + ckt, where ^ is the impact 
parameter of the light ray with respect to the origin of the coordinates, and fc is a unit vector directed along 
the direction of propagation of the light ray [32, section 7.3]. 

These equations are considered as fundamentals for creation of astronomical ephemerides of celestial bodies in the 
solar system [s^] and in any other localized system of self-gravitating bodies like a binary pulsar . In all practical 
cases they have to be extended to take into account the post-Newtonian corrections sometimes up to the 3-d post- 
Newtonian order of magnitude [gB*] . It is important to notice that in the Newtonian limit the coordinate time t of the 
gravitational equations of motion ([3]), Q coincides with the proper time of observer t that is practically measured 
with an atomic clock. In other words, the Newtonian limit of general relativity taken in the asymptotically-flat 
background spacetime equates the three time scales entering the Newton's law of gravity ([3]), the equation of light 
propagation (j4|) , and the equation of motion of an electron in atoms which solution defines the frequency of atomic 
transitions between different energy levels inside the atom that are used in practical realization of SI second [ssj . 

The post-Newtonian theory is well-established, mathematically elegant and self-consistent. It passed through 
numerous experimental tests with a flying color [65[- Is there any problem which we should worry about? The answer 
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is "yes" and the problem is that the background is not asymptotically-flat Minkowskian spacetime but is described 
by the FLRW metric, Qap. We live in the expanding universe. Therefore, the right thing would be to replace the 
post-Newtonian decomposition ([1} with the more adequate post-Friedmanian one 

gap = gap + ^ap , (5) 

where >iap is the metric perturbation around the cosmological background being represented as a series not only 
with respect to the post- Newtonian parameter 1/c but also, with respect to the Hubble parameter, 3-C. The latter 
terms appear in the post-Friedmanian series as a result of the Hubble expansion of the universe. Generalization of 
the post-Newtonian theory on the Minkowski background spacetime to that on the expanding background manifold 
is tantalizing and was a matter of considerable efforts of many researchers |1, [13, [13, [S, [3, [3, [s^ ■ A recent article 
[Toj summarizes the previous results and provides the reader with a number of other valuable resources. 

We notice that most of the previous works on celestial mechanics in cosmology used joining of two (for example, 
Schwarzschild and Friedmann) spherically-symmetric exact solutions of Einstein's equations, that was achieved in 



many different ways. McVittie's solution [45| is perhaps the most successful mathematically but yet lacks a clear 
physical interpretation [10]. What is still missing in literature is the precise mathematical formulation of the New- 
tonian limit for a self-gravitating localized astronomical system (a binary star, a solar system, a galaxy) having no 
specific spacetime symmetry while embedded to the expanding universe and coupled through the gravitational inter- 
action with the time-dependent background geometry. Theoretical description of the Newtonian limit for a localized 
astronomical system in expanding universe should conform with the results of the post-Newtonian theory obtained in 
the asymptotically-flat spacetime. Such a description will allow us to directly compare the equations of the "canon- 
ical" post-Newtonian celestial mechanics to its cosmological counterpart. Therefore, the task is to derive a set of 
the Newtonian-like equations in cosmology in some coordinates introduced on the background manifold, and to map 
them onto the set of the Newtonian equations ©-(HI) in asymptotically-flat spacetime. In addition, the equations of 
motion of electron in atom should be also re-formulated to take into account the Hubble expansion in order to check 
whether it appears in the equations and affects the accuracy and stability of the atomic time or not. All together 
such a theory of the post-Newtonian celestial mechanics and electrodynamics would be of a paramount importance 
for extending the tools of experimental gravitational physics to the field of cosmology and astrophysics, for example, 
to derive the cosmological extension of the PPN formalism [6^. The present article discusses the main ideas and 
principal results of such a theoretical approach. 

The problem of a paramount importance here is to find out whether the Newtonian limit of the cosmological 
perturbation approximates the Newtonian theory so that it contains only terms of the quadratic order in !K, or some 
terms being linear with respect to !K, survive. The same question should be answered for the orbital motion of 
electrons in atoms which, at least in principle, might be affected by terms being linearly proportional to the Hubble 
parameter. We will prove that the Newtonian theory for massive particles and the Coulomb approximation for electric 
charges can be indeed approximated up to the terms of the order of Jf^ in the same coordinate chart. However, this 
conclusion turns out to be inapplicable to the propagation of light. If equations of motion of electromagnetic signal 
are formulated in the same coordinates as those for massive particles and charges, they show the presence of terms 
which are linear with respect to 5f . The reason behind this intriguing difference is rather simple - Maxwell's equations 
for freely-propagating electromagnetic waves are conformally-invariant with respect to conformal transformation of 
the metric tensor while the equations of motion of massive particles and charges are not. Therefore, it is impossible 
to find out a single space and time coordinate transformation in FLRW manifold that will reduce equations of 
motion of particles, charges, and electromagnetic waves simultaneously to their classic form that they have in the 
asymptotically-fiat spacetime. We need at least two different time scales in order to do this. 

This simple mathematical fact leads to a fascinating theoretical conclusion - the Hubble expansion of the universe 
does not permit the existence of a single uniform time scale for physically-different gravitational and electromagnetic- 
wave processes. Time scale in gravitational equations of massive bodies known as ephemeris time (ET) [32|, coincides 
with the atomic time scale (AT) governed by the orbital motion of electrons in atoms giving access to SI unit of time 
[35} . Both time scales are identical with the cosmological time t of FLRW background universe. However, neither 
ET nor AT coincide with the time scale defined by a geodesic motion of a free electromagnetic wave under condition 
that the spatial coordinates defining positions of massive particles, charges, and photons belong to the same chart on 
the manifold. We show that the time scale of freely-propagating electromagnetic wave diverges quadratically as time 
passes on from AT and ET. 

The idea of using propagation of electromagnetic signal in space as a clock belongs to Marzke and Wheeler [i^ who 
proposed to run a light ray continuously between two mirrors which move along time-like geodesies of a background 
manifold. The electromagnetic time scale is measured by counting the number of reflections of light from the mirrors. 
In fact, the motion of mirrors may be not geodesic - it is sufficient to know their world lines with required precision 
to eliminate the effect of the Doppler shift on the frequency of light. Marzke- Wheeler's idea has been realized in 
frequency standards using a microwave cavity or an optical resonator. On the other hand, the same idea is used for 
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a long time in the Doppler-tracking radio technique but it was never associated with Marzke- Wheeler's optical clock. 
Doppler tracking of frequency of electromagnetic signal whose phase is locked to the primary atomic standard on 
Earth, allows a coherent measurement of the phase of electromagnetic wave propagating back and forth between a 
radio emitter on Earth and a spacecraft's transponder. Thus, the Doppler tracking is a high-quality radio resonator 
having the size of the solar system and formed by two "mirror" - radio emitter on Earth and transponder on board 
of the spacecraft. The above-mentioned quadratic divergence of the electromagnetic time scale from AT and ET is 
caused by the Hubble expansion, and the disagreement between the time scales accumulates proportionally to CKt^. 
If the divergence between the electromagnetic and AT/ET time scales is measured, it will allow to measure the local 
value of the Hubble constant !H without resorting to cosmological observations. 

We argue that the Pioneer effect has, in fact, a simple and natural explanation as the cosmological effect of quadratic 
divergence between the electromagnetic and atomic time scales governing the propagation of radio wave in the Doppler 
tracking system and the atomic clock on Earth respectively. It has nothing to do with the "anomalous acceleration" of 
spacecraft which was for a long time a matter of controversial theoretical speculations [37| . The present paper explains 
the Pioneer effect in the framework of general relativity without any resort to "new physics" or to an alternative theory 
of gravity. The most precise measurement of the local value of the Hubble constant can be achieved in a near future 
with a dedicated space mission equipped with precise atomic clocks and Doppler tracking system [s^ . The mission 
must be designed in such a way that minimizes various on-board generated systematic effects like heat, gas leakage, 
etc. (see section IXBj) 

In the following sections we describe the background FLRW manifold and its cosmological perturbation caused by 
a localized astronomical system. The perturbation is presented as the post-Friedmannian series with respect to the 
Hubble parameter IK. Only the linear terms will be found explicitly. The quadratic in 3-C terms are important for 
discussion of how dark matter and dark energy affects the orbital evolution of the localized system but they are not a 
primary concern of the present article and are given elsewhere (33| . The linearized cosmological perturbation of FLRW 
manifold is, then, decomposed into post- Newtonian series with respect to the parameter 1/c. The Newtonian limit of 
the perturbation is defined as an "osculating" solution which has the Newtonian form in the limit 1/c — > (X) but with 
the Hubble parameter being retained. The gravitational field obtained in this way, is used to derive the equations of 
motion of particles, and light. Maxwell's equations for electric charge distribution of an atomic nucleus are formulated 
on FLRW background with taking into account all terms depending of the Hubble parameter. Their solution is found 
in the form of the retarded integral which is expanded in series with respect to 1/c but with the Hubble parameter 
retained. Equations of motion of an electron orbiting the atomic nucleus are formulated in a covariant form and also 
decomposed in the series with respect to 1/c with the Hubble parameter included. We, finally, discuss the impact 
of the proposed theory on astronomical ephemerides and argue that it can explain the physical origin of the Pioneer 
anomaly 0, Q . The concluding remarks are about the local coordinates on cosmological manifold and the possible 
missions to measure the Hubble expansion in the solar system. 

II. NOTATIONS 

The notations used in the present paper are as follows: 

• Greek indices run through values 0, 1, 2, 3, and Latin indices take values 1, 2, 3, 

• Einstein's summation rule is used for repeated indices, 

• t and a;* = {x, y, z} are the coordinate time and spatial coordinates of the FLRW metric, 

• = {cri,x^} are the conformal coordinates with rj being a conformal time, 

• bold letters denote spatial vectors in Minkowski spacetime x = {x*}, 

• a prime /' = df/dt denotes a total derivative with respect to time t, 

• a dot / ~ df / drj denotes a total derivative with respect to the conformal time 77, 

• 9q = 9/9x" is a partial derivative with respect to a conformal coordinate x", 

• a vertical bar, denotes a covariant derivative of the function / with respect to the background metric, 

• a bar over function, /, means that the function belongs to the background spacetime manifold, 

• tensor indices of geometric objects on the background manifold are raised and lowered with the full metric ^q,^. 
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• tensor indices of geometric objects on the conformal spacetime are raised and lowered with the Minkowski metric 

fa/3, 

• the scale factor of FLRW metric is denoted a,s R — R{t), or as a = 0(77) — R[t{ri)], 

• the Hubble parameter J{ — R' /R, and the conformal Hubble parameter H = a/a, 

• the present paper operates with the WMAP value of the Hubble constant J{ = 71.0 ± 2.5 (km/sec)/Mpc [29l |. 
It corresponds to J{ = (2.3 ± 0.1) x lO^^^ s~^ 

The present paper deals with the theory taking into account all terms that are of the first order in the Hubble 
|)^ameter 'K. Terms which are proportional to Jf^ will be systematically neglected and published somewhere else 



III. BACKGROUND MANIFOLD 

The laws of nature evidently obey certain 
principles of symmetry, whose consequences we 
can work out and compare with observation. 

Steven Weinberg, Symmetry: A 'Key to Nature's 

Secrets 

The background manifold in the present paper is given by FLRW solution of Einstein's equations. The cosmological 
metric, gap, is spherically-symmetric and may admit the three-dimensional curvature of space which is characterized 
by a constant parameter k that take either of three values, k = {—1,0,1}. Cosmological observations consistently 
indicate that fc = [29| . and we accept this value in the present paper. The background cosmological metric with 
fc = is 

ds^ = -c^dt^ + R^{t) {dx^ + dy^ + dz^) , (6) 

where, c = 299792458 m-s~^, is the fundamental speed, R{t) is the dimensionless scale factor depending only on the 
cosmological time and = {x, y, z} are spatial coordinates. A continuous set of fiducial observers having fixed 
values of the space coordinates form the Hubble flow. Each Hubble observer measures the proper time t defined by 
dr^ = —c^ds^. For fixed spatial coordinates, t — t, which allows us to derive time t if the proper time r is measured 
with the help of some ideal clock. Usual assumption is that the proper time can be practically accessed with the help of 
atomic clocks. In expanding universe this assumption should be proved by deriving equation of motion of an electron 
in atom and demonstrating that the resulting equation can be reduced to the same form as in asymptotically-fiat 
spacetime. We provide such a proof in section rVIIIl 

It is also important to realize that the clocks of the Hubble observers are synchronized in the sense of the Einstein 
synchronization procedure which is based on the exchange of electromagnetic (light) signals [s^lg^]- Hypersurfaces 
of simultaneity for the Hubble observers are defined by the condition t = const.. The proper distance, £, on the 
hypersurface of simultaneity is £ — R{t)^Jx'^ + J/^ + z'^- 

For calculational purposes it is more convenient to introduce conformal coordinates x°' — {crj, x, y, z). The conformal 
time 77 is related to the cosmological time t by an ordinary differential equation 

dt , , 

where a{rj) = R[t{rj)] is the scale factor of the spatial part of the metric ^ expressed as a function of 77. The background 
metric gap in the conformal coordinates is conformally-flat. It may be worth mentioning that the conformally-fiat 
and asymptotically-flat spacetimes are not physically identical (50i] . 

gap = a^(?y)fa/3 ■ (8) 
FLRW spacetime manifold with the metric ([8]) has the Christoffel symbols, 

cf = - J{ {5"pu-, + 6^up - u^g^^) , (9) 

and the Ricci tensor 

C^Rap = ^' {gap - '^UaUp) + 3J{^5a/3 , (10) 
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where J{ = a' /a, = a~^(5o ^ four-velocity of the Hubble flow, and Ua — gapu^ — and a prime denotes a 

derivative with respect to the Newtonian time t. 

The dynamic evolution of the background manifold is determined by Einstein's equations 

SttG - 

Gap = —^Tap , (11) 

where 

Gap = Rap - l/2gapR% , (12) 

is the Einstein tensor, the Ricci scalar — g^'^Rf^i^, and 

Tap = (e + p) UaUb + pgap , (13) 

is the energy-momentum tensor of matter governing the dynamic evolution of the background universe. In FLRW 
universe this tensor has the form of the perfect fluid with the background values of energy density e, pressure p, and 
four- velocity Ua- The energy density is defined in terms of the rest mass density, p, and the internal energy per unit 
mass, n, as follows [47[ 

e = p(c'+n) , (14) 
where the internal energy density, 11, is related to pressure, p, by the law of conservation of energy (47| 

dll+pd = . (15) 

We assume that equations ([Tl]) -(fT5 |) can be solved for a particular equation of state and, hence, the time evolution 
of the scale factor a(r]) is well-defined. Since all geometric and thermodynamic quantities are the implicit functions 
of 0(77), their temporal evolution is also well-defined. 



IV. POST-NEWTONIAN APPROXIMATIONS IN COSMOLOGY 



We assume that a localized distribution of matter (a bare perturbation) characterized by the tensor of energy- 
momentum Tq,/3 is placed on the background cosmological manifold. It makes a localized N-body astronomical system 
where the bodies interact to each other through their own gravitational field. The presence of the localized matter 
distribution in addition to the background matter of the Hubble flow, perturbs the background geometry as shown in 
. Perturbed geometry is described by Einstein's field equations 

SttG 

Gap = {Tap + Tap) , (16) 

where Gap^ Tap, and Tap are the perturbed values of the Einstein tensor and the tensors of energy-momentum of 
the background matter and the localized system. The perturbations of the above tensors can be decomposed around 
their background values in a Taylor series with respect to the perturbation of the metric tensor >Cap, 

Gap = Gap + g''^^ + G|^^j + . . . , (17) 

Tap = fap+T^J^+Tj^^^ + ... (18) 
7ap = lap + ^2 + • • • , (19) 

where G^'j and T^^ are of the order of 0{>iap), cg] and T^J are of the order of 0(>f^^), and so on. The tensor Hap 

is the source of the linear perturbations and, hence, is the same order of magnitude as G^^'p and t'^^ . Perturbation 

back gravitational reaction of the perturbed background geometry on the structure of %ap ■ 
After substituting expansions ([5]), (fT7| -(fT9 l) to (fTH]) . and accounting for the background Einstein equations (fTTI) . 
we obtain a set of partial differential equations for the metric tensor perturbation Xap which iterative solutions yield 
a successive post-Newtonian (sometimes called post-Friedmannian [59|) approximations to the true geometry of the 
cosmological manifold Q. Notice that the background geometry is well-defined by the solution of the Einstein 
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equations (IT2t . We do not use any averaging procedure applied to (fT6)) in order to define the background manifold. 



This procedure is discussed, for example, in |3l| and may be important in higher-order post-Newtonian approximations 
but it has no impact on the results of the present paper. 

We focus on the Newtonian limit of the post-Newtonian approximations. The linear approximation is sufficient for 



this purpose. The field equations (I16|) in this approximation read 



ttG 



„2p(l) _ 



(t^S'+^-z^) , (20) 



where 

G 



''^P = - (lal3^^\^+ gal3B^\^~ Ba\l3 - Bfj\a^ , (21) 

;r complicated tensor ex 

perturbation 



is a rather complicated tensor expression depending on a particular cosmological equation of state [33| , the metric 



Lf} = -XaP + ^gaP^C^f, , (22) 

and vector ee P-^i^. 

It is well-known that an analytic coordinate transformation, 

x"" =x°'-C{x'^) , (23) 
induces a gauge transformation of the metric tensor perturbation (63 | 

I'a/S = f-a/S - Ca\l3 - ^I3\a + 9al3C\f ■ (24) 

Here, we keep only linear terms with respect to the gauge functions ^" assuming that these functions are of the same 
order of magnitude as ^c^p . Higher-order terms can be included in (j24p but this goes beyond the scope of the present 
paper (see [2^15^). 

Calculation of the transformed Einstein tensor perturbation, G'^^p , defined by equation pip along with perturbation 
([M)) . reveals that it is gauge-invariant in the sense that G'^^^ = G^^^. Perturbation, T^^, of the energy-momentum of 
the background matter is gauge- invariant, T^^^^ = , due to the background equations of motion of the background 



matter T°'^\p = [43,|54|. These gauge-invariant properties of tensors in the left side of pT|) demand that the tensor 
of energy-momentum of the localized system obeyed the covariant law of conservation, 'J"''|/3 = 0. This equation can 
be reformulated in terms of partial derivatives [64| 

dp + y=If "^^T^-^ = , (25) 

which is useful in defining the conserved mass, center of mass and a linear momentum of the localized astronomical 
system on the expanding FLRW manifold (see discussion in section IVip . We also notice that tensor T"*^ is gauge- 
invariant because it is the source of the linearized metric perturbation so that the gauge transformation (j24p changes 
its form in terms of the second order that are not considered in the present paper. 



V. POST-NEWTONIAN FIELD EQUATIONS IN COSMOLOGY 

We are going to solve the Einstein field equations (|20|) in the linearized Hubble-parameter approximation which 
means that all terms of the order IH^ will be discarded (these corrections can be found in [l^l). This approximation 
is sufficient for discussing the effects of cosmological expansion in the solar system. Notice that the linearized Hubble 
approximation also discards all terms which are proportional to the time derivative of the Hubble parameter, JC'. 
This is because the background Einstein equations (fTTj) suggests that W ~ Jf ^ [iO, . 

By inspection, we can prove that perturbation, t'^} — 0(3-C^), that is it has no terms being linear with respect to 
"K. Indeed, the tensor of energy momentum Tap of the background cosmological matter is that of the ideal fluid 



Tap = (e + p) UaUb + pgaP , 



(26) 
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which depends on the energy density e, pressure p, four- velocity of the matter it", and the metric tensor gajs- We 

expand tensor (^51) in Taylor series around its background value, Tap, as shown in ([TO]) . Perturbation T^^^ is found 
by taking a variation of Tap from all variables. It yields, 

'^aP = UaUpSe + {UaUp + (jap) Sp + {e + p) {UaSup + UpSUa) + pSgaP , (27) 

where e, p, u", and gap are the background values of the energy density, pressure, four- velocity and the metric 
tensor respectively; 5e, Sp, Sua and Sgap — Xap a-re the perturbations of these variables, and we assume the Einstein 
summation rule. The terms in (j27|) which are proportional to e and p, are of the order of in virtue of the 
background Einstein equations P^ - (IT^ . The perturbation 5e — {c^/cDpS^i, where Cs is the speed of sound of the 
background matter, p is the mass density of the background matter, and is the perturbation of the specific enthalpy, 
p, = {e+p)/p, of the background matter [13, El]- The e ~ pc^ ~ 0(J{^) due to the background Einstein equations. 
Finally, the second law of thermodynamics tells us that Sp = p6p ~ 0(!H^). These estimates allow us to conclude 
that --0{^{^), and can be dropped out from equation ([20]) for we consider only terms of the linear order in Di. 
It agrees with cosmological studies having been done by other researchers P, (see [3^ for more detail) . 
Equation becomes 

laP^'^lfj. + 9apB'^\^ - Ba\p - B p\a = ^4 Taff , (28) 

where ee I'^'^^p. It is gauge invariant and admits a gauge freedom (j24p in choosing the metric perturbations, lap, 
on the background manifold. The gauge freedom can be used in order to simplify solution of (l28l) . It is convenient to 
fix the gauge by imposing the following gauge condition [sj, [ssj 

oar 

= ^—l"Pup , (29) 
c 

which generalizes the harmonic (deDonder) gauge condition applied in the post-Newtonian approximations performed 
on the asymptotically-flat background [ij, [S^l , to cosmology. To understand the advantage of the gauge ([^ in the 
cosmological setting, we write down the covariant Laplace-Beltrami operator for the metric perturbation in terms of 
the partial derivatives as follows 

= \niap + —u^d^lap - - (^Ky'^U)]/? - - (J^^^^U)]. + O {%\ , (30) 

Qi C C C \ C / 

where 

^ 52 1 92 92 92 52 

□ = f ^ I I 

dxt^dx'^ c2 d?f dx^ d]p- dz^ ' 

is the wave operator in Minkowski spacetime, and we have neglected all terms of the order of IK^. Substituting ((29)) 
and (|30p to (|28p cancels out a significant number of terms, thus, reducing equation ()28p to a much simpler form 

2^ IGttG 
□/a/3 + —a^u^'d^.lap = ^j^a^T^^ . (32) 

This equation can be reduced to the wave equation by making use of the following steps. 
We define the conformal metric perturbations, ^ap and hap, by equations 

gap a^{r\){^ap^hap) , (33) 

lap = a'^ivhap , (34) 

so that 

haP = -lap + '^'iapl^fi , (35) 

and indices of hap and ^ap are raised and lowered with the metric fap- Substituting (p4l) to ([32l) and simplifying 
terms, yield 

2H IGttG 

□7a/3 2"'^';'^"'3 = —^aP , (36) 
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where the conformal Hubble parameter H = a /a = W /a and all terms of the order of have been discarded. 
Equation ([55]) explicitly contains a term depending on the Hubble parameter but it can be eliminated by introducing 
a new variable, a^ap- In terms of this variable, equation p6p becomes equivalent to 

□ (a^ats) = —^alap ■ (37) 

This is a wave equation in conformally-flat spacetime which can be easily solved with the standard technique of 
retarded potential if we impose the boundary condition of vanishing 07^^ and 07^^,^ at conformal past-null infinity 
[igj . The advanced potential violates causality at higher post-Newtonian approximations, and we do not consider it 
here (see [1^ for further discussion). The retarded solution of (|37p is 

AG f a(s)T„/3 is,x')d^x' 



aivhap{v,x) = ^ / f — '— , (38) 

Jv \x - x'\ 

where V is the volume of integration occupied by matter, Tq,^, of the localized astronomical system, and the retarded 
time s is defined on the null cone of Minkowski spacetime, 

s = i]--\x-x'\. (39) 

c 

Solution psp. p9p is given in conformal coordinates a;" = {crj, x"^}. It is Lorentz-invariant with respect to the Lorentz 
transformations of the conformal coordinates and satisfies the gauge condition (P^ . 

It is important to emphasize that the retarded metric potentials (j38p describe gravitational field everywhere both in 
the near zone and in the radiative zone of the localized astronomical system emitting gravitational waves. Propagation 
of gravitational waves in the expanding universe takes place on the hypersurface of a null cone (1391) in the conformal 
spacetime. Therefore, the speed of weak gravitational waves is equal to the speed of light, c, referred to the conformal 
coordinates a;" = {cr;, a;}. 

Solution of (l37l) is defined up to the residual gauge transformation ((23)) with the gauge functions obeying the 
homogeneous wave equation 

a{aC) = . (40) 

This equation is valid up to the terms of the linear order of 0{3{). Quadratic corrections of the order of O(Jf^) to 
equation ([40]) are given in [ssj . 



VI. THE NEWTONIAN LIMIT OF THE FIELD EQUATIONS 

Equation (|38[) yields the post-Newtonian potentials jap in the linearized approximation with respect to the Hubble 
parameter J{ which is incorporated implicitly to the solution through the time-dependent conformal factor a = a{vi). 
Only the potential 700 is relevant in discussion of the Newtonian-theory approximation which is defined as a limit of 
the retarded solution (p8l) when the parameter 1/c — >■ 0. 

The Newtonian potential, U , in expanding universe is defined as C/ = — c'^7oo/4. In what follows, we shall focus 
on the localized system consisting of one massive body. Our treatment can be easily extended to N-body problem 
for, in the approximation under consideration, the Newtonian potential obeys the principle of a linear superposition. 
Equation ([38)) yields 



G r a(s)Too is,x')d^x' 
Jv a{ri)\x-x'\ 



uir,, x)=-j , (4i) 



where the integration is performed over the volume of the massive body. We assume the matter inside the body moves 
slowly, ^ c. Hence, we can expand the retarded function in the integrand of ()41[) in the Taylor series around the 
present time 77. It leads to 

, G f ToQiv,x')d^x' G f d , , , ,,,3 ^ f G\ 



The first term in right side of (jH)) looks exactly as the Newtonian potential of a conventional Newtonian theory. Yet, 
we cannot use this expression for doing the multipolar expansion with respect to the gravitational harmonics - we still 
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need to define the concept of the conserved mass, dipole moment and other multipoles of the astronomical system 
that will be independent of the Hubble parameter. Moreover, the correct Newtonian limit demands that the second 
integral in the right side of (H^ would be, at least, of the second order of magnitude with respect to M which we 
systematically discard. 

These issues can be resolved with the help of the covariant law of conservation of the tensor of energy-momentum 
'^ai3 - equation (f25|) . Making use of ([9]) for the ChristofFel symbols results in a more explicit form of (|25|) . 



19, , d . H 

-— (aXoo) - TT- (a^oj) = aTkk , (43) 

c or] ox3 c 



Id U 

-—{a%o)- — ia1^j) = aT,;o , (44) 

c or] 0x3 c 



Implementing (|43| for calculating the second integral in (|42|) . yields 

/ |- (aToo) d^'x = -aH f ^kkd'x , (45) 
Jv or] Jv 

where we have applied the Gauss theorem to eliminate the integral for the divergence of the vector field aToj • Further 
application of the laws of conservations P5)) . pij) brings about the virial relationships 



Jv 2c drj 



(46) 



a%od^x = / aTooa;^^^ + 0(5{) , (47) 

v c drj 

where 0{K) indicates the residual terms of the order of the Hubble parameter . 

Let us assume that distribution and motion of matter inside the localized system, is stationary. Relationships (j46l) , 
(l47l) suggest, then, that the right side of equation p5)) is of the order of and can be neglected in the linearized 
Hubble approximation. This result also suggests that we can define the conserved mass, M, of the massive body 
residing in the expanding universe by the following integral 

M = 1 / air])Too{v,x)d^x . (48) 
c Jv 

The mass is conserved, dM/drj — 0, up to the terms of the second order in Jf^. Since d/dr] = a{rj)d/dt, the law of 
conservation of mass can be recast to 

where we have used the cosmological time t that coincides with the proper time measured by the Hubble observer. 
We define the dipole moment, P , and the linear momentum, P*, of the massive body in the expanding universe by 

r ^ \ f a{r])%oo{v,^yd^x , (50) 
c Jv 

P' ^ f a^{r])%o{r],x)d^x . (51) 
c Jv 

Equation (|T7)) yields a relationship between a time derivative of /* and P*, 

dP 

Taking a time derivative from and applying (1441) along with (j47p yields the law of conservation of a linear momentum 
of the isolated self-gravitating body in the expanding universe 

dP' , . 

Equations (|52p , (15^ tell us that we can always chose the inertial coordinates in the expanding universe such that the 
center of mass of the body is defined by the condition, P = 0, and it remains at rest or moves with a constant velocity 
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in these coordinates. After choosing the center of mass at the origin of the inertial coordinates, a;" = {crjjX^}, the 
muhipolar expansion of the Newtonian potential (|42p reads 

C/(77,a^) = 4^^ + ... , (54) 

where eUipses denote quadrupole, octupole, and higher-order harmonics, and |a;| is the radial distance from the spatial 
origin of the conformal coordinates to the field point. The multipolar expansion (j54[) is valid in the linearized Hubble 
approximation with all terms of the order of having been abandoned. 

The Newtonian potential (|54|) looks like the potential of a point-like mass M with the effective gravitational constant 
25(77) = G/a(r]), depending on time. It may lead to a conclusion that the expansion of the universe encompasses a 
secular decrease of the gravitational constant - the idea pioneered by P.M. Dirac [16] and reiterated by R. Dicke in 
his weak anthropic principle [l3] and a scalar-tensor theory of gravity (along with C. Brans) 9]. This conclusion is 
too straightforward and is invalid in general relativity. The reason is that the gravitational constant can be measured 
only through observing orbital motion of celestial bodies and processing the observational data on the basis of the 
Newtonian equations of motion. Scrutiny analysis given in next section shows that the scale factor a{r]) drops out of 
the equations of motion after making a transformation from the conformal to local coordinates, so that the effective 
gravitational constant 25(77) gets rescaled and its variation cannot be observed - only the constant value of G can be 
measured. 



VII. THE NEWTONIAN EQUATIONS OF MOTION IN EXPANDING UNIVERSE 

Astronomical observations make use of test particles of two sorts: with a non-zero rest mass and massless photons 
(light, radio). In order to understand the impact of cosmological expansion on the motion of particles we have to 
derive relativistic equations of their motion and to analyze their Newtonian limit which is reached for the value of 
the parameter 1/c — > 0. The task is to find out the analytic dependence of the Newtonian equations of motion on the 
Hubble parameter 'K. 

Any test particle, if it is unaffected by an external force of non-gravitational origin, moves along a geodesic worldline 
which equation is given by 

^ + rv^^^^o. (55) 

where r is an affine parameter on the world line of the particle, T"^^ is the perturbed affine connection given by 

r"pi. = ^5"^ (5/3M,'' + 5/3.^,^ - dt^i^.fi) , (56) 

and the perturbed metric gap is defined by p3p with taken from psp and (j38p . The affine parameters for time- like 
and light-like geodesies are different _64J. Irrespectively of the nature of the affine parameter we can always choose 
the coordinate time 77 as an independent parameter along geodesies. Equation (|55p parameterized by the coordinate 
time 7;, reads 

+ rV.i'^i" - ^T\,x''x''x' = . (57) 

We elaborate on these equations in the next two subsections. We shall focus primarily on discussion of the Newtonian 
limit of these equations in application to test particles with a rest mass, and to light (massless particle, photon). 



A. Equations of motion of test particles with non-zero rest mass 



Calculating the Christoffel symbols and substituting them to ([57)1 , yield the equations of motion of test particles 

GA4 

X = -Hx ; — -ttX , (58) 

where all terms of the post- Newtonian order ^ 1/c^ and the order of , have been omitted. These equations 
are written in the conformal coordinates {77, a;} and apparently have terms which makes them looking unlike the 
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Newtonian equations of classical mechanics. In particular, we notice the presence of the conformal Hubble parameter 
H and the time-dependent scale factor a{'q) that again makes an impression that the effective gravitational constant 
25(ry) = G /a{'q) appears in the equations. 

It is remarkable that there are transformations of the conformal coordinates which convert (j58p to exactly Newtonian 
form of the classic mechanics. The coordinate transformations are 



t = i a{r]')dr]' , r = a{-q)x , (59) 

where we integrate from the initial epoch, rjQ , and assume that at the initial epoch the dimensionless conformal factor 
o-ivo) — 1- These transformations reduce equation (|58p to the exact form of the classic Newtonian equations of motion 

GM , , 

r" = -—r , (60) 

where r — |r|, and a prime denotes a derivative with respect to time t. 

The obtained equation is fully identical with the Newtonian equations of motion in the classic celestial mechanics, 
and does not predict any "cosmological force" that would be proportional to the Hubble constant 'K. All terms 
depending on H canceled out after applying the conformal space and time transformation (|59p . The explicit de- 
pendence of the Newtonian force on the scale factor, 0(77), has disappeared as well. Our derivation of equation ([SO]) 
helps to understand the result obtained by Brumberg and Krasinsky with a different mathematical technique. 
They analyzed circular equations of motion by solving equations (j58p explicitly, and arrived at the conclusion that 
the circular orbits are not affected by the cosmological expansion. In fact, this is a consequence of the conformal 
transformation (j59p which allows us to prove that in FLRW universe any massive particle moves in accordance with 
the Newtonian equations of motion without any presence of the Hubble parameter H. More advanced post-Newtonian 
analysis of the residual terms in the Newtonian equations of motion show that the residual terms are proportional 
to ~ Hx^/c and ~ H^. Hence, any Keplerian orbit of massive particle (not only circular) remains unaffected by the 
cosmological expansion. This conclusion has been validated by other researchers who used different mathematical 
techniques [j. [lOl, lii] . 



B. Equations of motion of light 

Light is a wave of electromagnetic field which obeys homogeneous Maxwell's equations in free space. These equations 
are conformally-invariant as we prove in the next section. For these reasons, the equation of motion of light in the 
conformal coordinates {77, x*} is the same as in Minkowskian spacetime (remember that we consider the case of FLRW 
universe with space curvature fc = 0) 

x = 0. (61) 

Hereby, we have neglected the gravitational force exerted on light by the localized astronomical system as it produces 
the deflection of light and time delay which are of the post-Newtonian order and are not of the primary concern in 
the present paper. We draw attention of the reader that equation (|6ip does not have any residual term depending on 
the Hubble parameter IH. This is because of the conformal invariance of the Maxwell equations [g^l . 

Equation (ISTt can be also derived from ([57]) by calculating the Christoffel symbols, making use of equation of light 
cone, f = 0, and neglecting the gravitational deflection of light caused by the localized astronomical system. 

Equation (|6ip tells us that light in conformal coordinates of the conformally-flat FLRW manifold propagates with the 
constant speed c along straight lines. Equation (|5T|) looks exactly the same as the Newtonian equation of propagation 
of light but one has to use this analogy with care. The point is that the conformal time, 77, and space coordinates, 
X, entering (|6T|). are not those {t, r} used in the Newtonian equations of motion ((60)) of massive particles. Therefore, 
we can not use solution (1611) to analyze the orbital motion of massive particles. 

In order to make equations (j61l) compatible with the Newtonian equations of motion (|60p we should make use of 
the conformal coordinate transformations (|59l) in equation ([CT|) . It recasts (|6ip to 

r" = ■Xr' . (62) 

This equation should be compared with ^ that we would expect in the Newtonian limit of general relativity in 
coordinates, {t,r}. However, the conformal nature of the cosmological spacetime and the conformal invariance of 
Maxwell's equations does not support this expectation and bring about equation (1621) that yields the correct Newtonian 
limit for the equation of propagation of light in cosmology with the space curvature fc = 0. Similarly to (|4]), equation 
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((62)) describes propagation of photons along straight hnes but we notice the presence of a force, ?fr', in its right side 
that is acting on photon due to the expansion of the universe. Appearance of this force points out that coordinates, 
r}, can not be considered as inertial coordinates for description of hght propagation. 

Equation (p^ can be solved exactly but it will introduce terms of the order of and higher, which we neglect. 
We are interested only in keeping linear terms with respect to IK. For this reason, we can assume that J{ in the right 
side of (|62l) . is constant. It allows us to directly integrate this equation with respect to time. The integration yields 
a coordinate velocity of light, 

r' = (1 + m)ck , (63) 

where the unit vector k = {fc'} is directed from the point of emission of light toward the point of its reception 
{Sijk'^y = 1). The non-inertial force 3-Cr' accelerates photons with respect to spatial coordinates, r, so that the 
coordinate speed of light, r', does not remain c but grows directly proportional to the cosmological time t. We 
emphasize that there is no contradiction with special relativity because the physical speed of light, measured locally 
at time t = 0, remains equal to the fundamental speed c. One more integration of ([63)1 yields trajectory of the light 
ray in coordinates {t,r}, 

r^$ + ck(^t + ^^t^^ , (64) 

where ^ = {^*} is the (constant) impact parameter of the light ray with respect to the origin of the spatial coordinates. 
The impact parameter is orthogonal to k in the sense that the Euclidean dot product of the two vectors vanishes, 
Sijk^S,^ = (more details about the geometric properties of the parametrization see in our textbook [s^, section 

Another way to understand the propagation of light in the expanding universe is to apply a slightly different 
conformal transformation of time and space coordinates, namely, 

^ = a{v){v - Vo) , r^a{'q)x, (65) 

where A is a new time coordinate which differs from the conformal time rj by the time-dependent scale factor a{r]). 
This transformation leaves equation (pT|) form-invariant, in the sense that if all terms of the order of Jf^ are neglected, 
the equation of light propagation in coordinates {A, r} has the same form as in the conformal coordinates {rj, x}, 



dX 



2 . (66) 



Equation ()66|) means that the coordinates {A, r} can be considered as inertial coordinates for describing propagation 
of light. 

Solution of equation ()66|) is a straight line 

r = I + cfcA , (67) 

with the constant coordinate speed of light, dr/dX = ck, and where k and ^ have the same meaning as in (|64p. By 
expanding the conformal factor in (1651) in a Taylor series around rj — rjo and making use of (|59p , we get a relationship 
between the two time scales, 

X = t + i , (68) 

that matches equation (|57| with equation (|64p . 

The lesson we have learned from the analysis of the equations of motion in the conformal FLRW spacetime, is 
that there is no single inertial coordinate chart in the expanding universe in the Newtonian description of motion 
of massive particles and light. Transformation to inertial coordinates, {t,r}, eliminate all Hubble-dependent terms 
of the order of !K from the Newtonian equations of motion of massive particles but it can not eliminate such terms 
from the Newtonian equation of motion of light. On the other hand, transformation to the inertial coordinates {A, r} 
eliminate all Hubble-dependent terms of the order of Jf from the equation of motion of light but put them back to 
the equations of motion of massive particles. 

There is a physical difference between the two time scales, t and A. The cosmological time t coincides with the 
proper time of the Hubble observers while the time scale A diverges quadratically from t and does not coincide with 
the proper time. Special relativity postulates that proper time can be measured with an ideal atomic clock. This 
postulate should be checked in cosmological spacetime to make sure that atomic clocks do measure the proper time. 
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VIII. RELATIONSHIP BETWEEN ATOMIC AND PROPER TIME IN COSMOLOGY 



General relativity postulates that the proper time of observer is measured by an ideal clock carried by the observer. 
The most precise practical realization of the ideal clock is the atomic frequency standard which operational principle 
is based on the use of electron's energy level transitions in atom to produce a stable frequency that drives the clock 
(4]| . In classic regime the atomic frequency is defined by the difference between the energies of two orbits of an 
electron moving around nucleus of the atom in accordance with the second Newton's law with the central Coulomb 
force. Additional Bohr's postulate is required in order to match the classic motion of electron in atom with quantum 
mechanics but it is not essential for further consideration. Thus, in order to make sure that the proper and atomic 
time scales in the expanding universe coincide, we have to prove that the orbital motion of an electron in a hydrogen 
atom is not affected by the expansion of the universe, at least, in the linearized Hubble approximation. The problem 
was studied in 0, [13, EH and the authors have reached a consensus that the expansion of the universe should not affect 
the stability of atomic frequencies. Here, we refine the proof of this statement with a more rigorous mathematical 
technique. 

Electromagnetic field tensor Fa/3 = daAjs — dpA^ is defined in terms of a vector potential Aa- General-relativistic 



Maxwell's equation for electromagnetic vector A" in a curved FLRW spacetime are [47 1 



— A-TT 

Aj^\^- A^\f3a- Ro.fi A^ = -— Jo., (69) 

where Ja is a four- vector of electric current. Usually, the Lorentz gauge condition, A^^p = 0, is imposed on vector- 
potential in order to eliminate the middle (gauge-dependent) term in the left side of equation (j69p and reduce it to 
the de Rham equation j47]. Unfortunately, such a gauge condition does not simplify the last term in the left side of 
()69p depending on the Ricci tensor. We have discovered that another gauge condition is more elegant (compare with 
([29| for gravitational case), 

A'5,^ = -^A%. (70) 

This gauge condition has been also found by M. Ibison [l^, in his study of electrodynamics in cosmology with a 
future conformal horizon. 

We substitute ([70]) into the second term and use expression (fTU]) for the Ricci tensor. It turns out that after 
performing the covariant derivatives and reducing similar terms, all the terms depending on the Hubble parameter 
cancel out in all orders, and we arrive to exact Maxwell's equation in FLRW spacetime (with space curvature fc = 0), 

DAa = aV„ , (71) 

c 

where □ is the wave operator in the Minkowski spacetime defined in (PT|) . It is remarkable that it look similar to 
special relativity in Minkowski space-time except for the scale factor, a = a(r/), that appears in the right side of (|71l) . 
The reason for the simplicity of (|7ip is the conformal invariance of electromagnetic field which makes the left side 
of (|69|) be equal to the left side of ((7T|) up to the conformal factor 1/a^ 63]. We take a physically-relevant retarded 
solution of equation ([71]) . 

where the retarded time s is defined in p9l) , and the volume integral is performed over the charge distribution inside 
the atomic nucleus. Because the retarded argument s in the solutions for electromagnetic and gravitational field is 
defined by the same equation p9|) , weak electromagnetic and gravitational waves propagate with the same speed c in 
the conformally-flat space of FLRW universe. 

We are interested in the Coulomb approximation for the potential Aa to describe the electric field produced by 
atom's nucleus. In this approximation, it is sufficient to consider only the electric potential (j) = cAq. Everywhere 
inside the near-zone of the atom, the retarded argument of the vector potential can be expanded in a Taylor series 
around time 77. For the electric potential it yields 

f a^ir,)Mr,,x')d^x' , 1 ^ 9 r 2, ^ w '^]^3>^nf^\ (7'^\ 
^('^' ^) = Jy \^\ +cjyd^^'' ^"^^-^'^"^^ ^)]d-+0[^-^y (73) 



Taking covariant divergence from both sides of equation (|69l) reveals that the divergence of the left side vanishes 
identically. It means that the electric current is conserved, 

= , (74) 
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which is equivalent to 



The total electric charge Q is defined by 



Q = - a^{7])Joi7],x)d^x . (76) 
Jv 

The charge is conserved, dQ/di] = 0, up to terms of the order of Jf^ as follows from ([75)) . Replacing the derivative 
with respect to the conformal time 77 by the derivative with respect to cosmological time, the law of conservation of 
charge in expanding universe becomes 

dQ/dt = 0. (77) 

Because of charge's conservation, the second term in right side of (j73p vanishes while the first term can be expanded 
in terms of the electric multipoles - the constant charge Q, the dipole electric moment Qi, and so on, 

In what follows, we neglect the dipole and higher-order multipole moments of electric field. 

The atomic frequency of electromagnetic wave emitted by atom, is defined by the difference between two values of 
orbital energies of an electron when it jumps from one orbit of the atom to another. The orbital motion of an electron 
in a curved spacetime is given by the second Newton's law with an electromagnetic Lorentz force in its right side [131 
while the left side of this law is a covariant derivative from the linear momentum of electron = mdx°' /dr, 

d^x- dx^dx^\ dx^ 

'-^1^)='^ ^^' (^^^ 

where r is the affine parameter along the electron's orbit, e is the charge of electron (e < 0), and m is electron's 
mass. We accept that both charge and mass of electron remain constant in expanding universe in accordance with 
the definitions of these conserved quantities given in (|76p and The Christoffel symbols in (175)) account for the 

presence of gravitational field. It is sufficient to take into account only the background value of the Christoffel symbols, 
r"/ii^, as the metric tensor perturbations cause negligible effects at least, in the solar system. Strong gravitational 
field of a localized gravitational system (a black hole) may affect the clock's performance 21] but we do not discuss 
the strong field effects in the present paper. 

If the conformal time, 77, is used for parametrization of electron's orbit, equation of motion (j79p assumes the following 
form 

+ f V^i^i'' - if ° ^i^'i-'i' ^ — (f'b- -F%x'] x^t . (80) 
c m \ c J 

In the slow-motion approximation f = 0(77), F^q = Fio/a^, and the electric field of the atomic nucleus is Ei = cFio = 
di4>. By neglecting relativistic corrections of the order of 1/c^, we obtain the equation of motion of electron in the 
expanding universe 

x = -Hx + ^^, (81) 
ma \x\'^ 

which mathematical structure is similar to the Newtonian equation of motion of massive test particles (|58l) . Hence, 
by doing conformal transformation ([59]) we can reduce ([8T|) to the classic form of the second Newton's law with 
Coulomb's force 

mr" = . (82) 

Thus, choosing coordinates {t, r} on the expanding cosmological background eliminates any dependence of the orbital 
motion of electrons on the scale factor and the Hubble parameter. This proves that the atomic frequencies are not 
affected by the expansion of universe, at least, up to terms of the order of 3^^. Furthermore, the time parametrization 
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of equation (I82p makes it evident that the atomic time (AT), measured by the ideal atomic clock carrying out by a 
Hubble observer, is identical with the proper time, r, of the observer as well as with the cosmological time t. 

Comparing two equations, (|60p and (|82p . we observe that the parametrization of the orbital motion of massive 
particles in gravitational field and that of electrons in atoms is done with the same parameter - the cosmological 
time t. The argument of the gravitational equations of motion of massive particles ((60)) is called the ephemeris time 
(ET). Our consideration makes it clear that in general relativity ET can be measured in two different ways: (1) by 
observing motion of massive particles in gravitational field, and (2) making use of atomic clocks. If general relativity 
is valid we should not see any divergence between the two time scales - ET and AT. 

In contrast to ET and AT time scales, the time scale. A, that parametrizes the equations of motion of photons 
in expanding universe, diverges quadratically from the atomic time (AT) in accordance with equation (1681) . This 
divergence has important physical consequences in astrophysical applications and in testing general relativity as it 
allows us to measure the local value of the Hubble parameter J{ independently of cosmological observations. 



IX. THE IMPACT ON ASTRONOMICAL EPHEMERIDES 



Modern ephemerides represent possibly the 
greatest dynamical system of all time... 

Myles. E. Standish, [Hi] 



A. What equations should we utilize for space navigation? 

Astronomical observatories and the national Space Navigation Centers (like NASA JPL) operate with the computer 
codes calculating astronomical ephemerides of celestial bodies under assumption that the Newtonian equations ([2])- ([4]) 
are physically-adequate and applicable on equal footing to the equations of massive particles and light in the inertial 
coordinates {t, r} whose origin is chosen at the barycenter of the solar system. Post-Newtonian corrections to the 
Newtonian equations are included to the computer codes to reach the computational precision being comparable 
with the current accuracy of astronomical observations [i^. The post-Newtonian corrections in the solar system 
ephemerides have been thoroughly studied and are well-known [32] . Hence, we do not discuss them in the present 
paper. We shall also neglect the motion of the solar system with respect to the cosmological reference frame defined 
by the congruence of the Hubble observers, and assume that the worldline of the solar system barycenter has fixed 
cosmological coordinates. It allows us to identify the cosmological time, t, with the barycentric coordinate time (TCB) 
used in the ephemeris astronomy [57.] . Thus, the hypersurface of simultaneity of the cosmological time t coincides 
with that of TCB when we use the spatial coordinates r — a{rf)x. 

Fitting astronomical observations to theoretical predictions utilizes the set of the Newtonian equations which leading 
terms are 

GM , , 

r" = -—r , (83) 

= r2i . (84) 

Here the first equation (|83|) is used for calculating coordinates, r = r{t), of planets and spacecrafts orbiting the Sun 
with Af being the mass of the Sun. The second equation (|84p is the light travel time (light-cone) equation used 
for relating the radar distance, £21 = c{t2 — ii), to the coordinate distance, r2i — \r2 ~ ri\, between the point of 
emission of light with coordinates {ii,Ti}, and the point of reception of light at the point with coordinates {^21 ''"2} • 
Light-cone equation (|84p can be viewed as providing a set of initial and boundary conditions for equation (j83p . In 
principle, a rigorous definition of the radar distance is based on the time interval a light signal takes in order to travel 
from observer to a celestial target, and back. The light round-trip time is a proper time, r, measured by observer 
equipped with atomic clock. Therefore, strictly speaking, in addition to equations ()83p . (|84p . we have to write down 
one more equation that transforms the proper time of observer, r, to TCB—t. This equation is r = t + 0(c~^), and 
is called the time ephemeris [3^ . Its precise formulation and solution with the account for relativistic corrections, 
^ c~^, are well-known and has been discussed, for example, by Fukushima Relativistic corrections in the time 
ephemeris equation are not important for discussing the Newtonian limit in cosmology. Hence, we omit them in the 
time ephemeris equation and equates t = t, which is interpreted as the proper time of the Hubble observers. 

The inertial coordinates {t,r} used in (|83| . (l84t have direct physical meaning. Indeed, time t =TCB enters 
gravitational equations of motion (|83p and defines the ephemeris time ET which can be compared with the atomic 
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time (AT) measured by observer on Earth. The distance r2i is known in cosmology as the proper distance between 
two points lying on a hypersurface of a constant time t If the two points are infinitesimaUy close to each other, 
the proper distance r2i coincides with the radar distance £21 [38| in accordance with equation (j84p which is accepted 
as exact in the process of calculation of ephemerides. However, this equation is not true in the expanding FLRW 
universe. The distances r2i and £21 differ from each other, if the two points are sufficiently far away and measuring 
the radar distance requires integrating the light-ray equation (|62|) along the light cone. Indeed, the Minkowskian 
structure of the light cone in expanding universe is preserved in the light-ray inertial coordinates {A, r}, but not in 
{t,r}. Therefore, relationship (1841) between the radar and proper distances is only an approximation which diverges 
quadratically with increasing the light travel time between the two points. More specifically, equation (|84|) should be 
replaced with 

c(A2- Ai) =r2i , (85) 

where 

Ai = ti + ^'Ktj , X2 = t2 + ■ (86) 

Time scale A, is uniform for description of propagation of light while the time t is a uniform time scale for celestial 
dynamics of massive particles. Current practice ignores the difference between the time scales, A and t. It is important 
to evaluate the consequences of this ignorance as it can lead to a secular divergence between predicted and observed 
positions and velocities of celestial bodies. It can be misinterpreted as an gravitational effect of a yet unknown force 
in the equations of motion of celestial bodies leading to "anomalous orbital acceleration" 0, H, [s^l • We shall analyze 
the consequences of replacement of an approximate equation (j84p for exact equation (j85p in case of two astronomical 
techniques - radar/laser ranging and the Dopplcr tracking. 



B. Range Measurements 

We use equation ([BS]) and express times ti, t2 in terms of the coordinates ri, r2 of the points of emission and 
reception of the light ray by making use of equation Taking the Euclidean dot product of the vector equation 

(|64p with itself and accounting for the orthogonality of the impact parameter ^ and the unit vector fc, yield for the 
radar distance £21 = c{t2 — ti) the following result, 

£21 = r2i + 5r2i , (87) 

where 




Equation ()87|) makes range measurement more precise, as compared with equation (|84[) . by introducing correction, 
(5r2i, to the proper distance, £21, due to the expansion of the universe. This correction is fairly small and amounts 
to (5r2i ~ 10~^ cm in case of the lunar laser ranging, and W21 — 0.1 cm in case of Earth-Mars radar ranging. These 
numbers elucidate that the expansion of the universe would not call the shots in the astronomical ephemerides if 
only the radar measurements were used in fitting theory to observations. However, the astronomical ephemerides 
also utilize another technique known as the Doppler tracking. This technique measures velocity of spacecraft which 
is equipped with a high-stable clock and a transponder allowing us to continuously track the phase of the Doppler 
signal for sufficiently long interval of time. 



C. Doppler tracking 

Doppler tracking is a method of measuring and control of deep space probe's velocity by making use of a radio- 
communication system onboard of the spacecraft via radio transmissions from Earth to spacecraft (the uplink) and 
returning telemetry to the ground via transmission from the spacecraft to the Earth (the downlink) |4^ . Analogue 
of the Doppler tracking in cosmological observations is the measurement of the frequency shift of a spectral line of 
light emitted by atoms of distant objects - galaxies and quasars - that led Edwin Hubble to the discovery of the 
cosmological expansion of the universe. The Hubble expansion causes shifting of the spectral lines of the light from 
receding cosmological objects toward the red part of electromagnetic spectrum (decrease in frequency) as compared 
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with the laboratory spectrum. Measuring the amount of the red shift astronomers can determine a radial velocity, v, 
of the object by making use of equation for the Doppler shift [13, [IJl. Cosmological distance, r, to the object relates 
to its radial velocity v by the Hubble law which, in the first approximation with respect to the Hubble parameter, is 

v = 3<r . (89) 



The distance r entering equation (|89p is the proper distance similar to definition of r2i in the range equation (j84p . 
There are other definitions of distance in cosmology based on photometric, parallactic or other kind of measurements 
[13] but they coincide in the linearized Hubble law approximation (|89|). Spacecraft's Doppler tracking in the solar 
system uses a coherent measurement of phase of electromagnetic wave transmitted from Earth to the spacecraft and 
re-emitted back to Earth. The phase and frequency of the electromagnetic wave - both transmitted and received back 
- are continuously recorded by the Deep Space Network (DSN) tracking system. The Doppler tracking represents an 
ideal radio-communication system for performing the most precise celestial mechanical experiments with spacecrafts 

We can comprehend how the cosmological expansion affects the Doppler tracking measurements by calculating the 
Doppler effect with the help of equation (|85l) for the light travel time between spacecraft and observer on Earth. We 
denote the proper time of the atomic clock on the DSN tracking system as ti, and the proper time of clock on-board 
of spacecraft as T2. The proper time t as a function of the conformal time rj, conformal coordinates x and conformal 
velocity x°- = (c, x), of observer/spacecraft, is defined by the integral relationship 

y'l - c-2 (a;2 + hapx°'xf')a{ri')dr]' , (90) 

that follows directly from the perturbed cosmological metric ([55]) and definition of the proper time, c^dr^ = 
—Qafii^'i^drf'- We shall neglect the special relativistic effects and those due to the gravitational perturbation hap 
which are proportional to l/c^. These effects are important in high-precision space navigation but are inessential in 
our discussion of the influence of the cosmological expansion on the frequency of a radio wave bouncing up and down 
in the radio-communication link connecting Earth and spacecraft. In this approximation the proper time ri of the 
terrestrial observer and that of spacecraft, T2, coincide with the TCB times of emission, ti, and reception, ^2 , of the 
signal. The post- Newtonian corrections in the transformation between r and t are discussed in detail in [32l |. 

Let us consider, first, the up-link Doppler frequency relationship, when a radio wave is transmitted to spacecraft 
from the ground-based DSN antenna. The frequency of the transmitted radio wave, vi = 1/dri, where dri is the 
period of the wave. The frequency of this wave, when it reaches the spacecraft, is z^2 — ^/dT2, where dT2 is the period 
of the received wave. In the approximation when all terms of the order l/c^ are neglected, we have dri = dti and 
dT2 — dt2. The transmitted and received frequencies are linked by the light-cone equation ([55)) . Taking the total 
differential of this equation with respect to dti and dt2 and accounting for (j86p , yields 



_ dh _ l + 0<t2- fc21 ■ f32 
dt^~ l + mi- fc2i • /3i 



(91) 



where the up-link unit vector k2i — T2i/'"21 is directed from Earth to the spacecraft, l3i — Vi/c, (32 — i'2/c, 
Vi = dri /dti is a velocity of the tracking station at time ii, and V2 = dr2/dt2 is the spacecraft's velocity at time 
t2. These velocities are modelled by solving equations of motion (|55)) with the appropriate initial conditions jisj . 
Equation (lOT)) can be revamped by making use of the property of the unit vector: fc2i • ^21 = 1. It yields 

^ 1 - fc21 ■ {(32 - fc2lTO2) 
1^1 l-fc21-(/3l-fe2lJ{ti) ' ^ ' 

which shows that the Hubble expansion modifies velocities of observer on Earth, ^i, and spacecraft, ^2 by subtracting 
a small (blue shift) correction IKt directed along the line-of-sight of observer toward spacecraft. 

Let us assume the Earth-spacecraft distance is about the size of Pluto's orbit that is ~ 40 astronomical units. Then, 
the light travel time <2 — — 2 x lO'' s, and the Hubble expansion yields a fractional uncertainty in the received 
frequency of the order Si^/i^ = ^{t2 — ti) ~ 5 x lO"^"*^ which is too small for detection in a single Doppler frequency 
measurement. However, the Doppler tracking is not a single measurement of frequency. It is actually a continuous 
series of coherent measurements which consists of a consecutive exchange by a radio wave bouncing between the 
DSN transmitter and a transponder on-board of spacecraft. The radio-tracking system observer-spacecraft makes up 
a radio-resonator with a high-quality factor, Q, and a radio wave being phase-locked inside it. Thus, the Doppler 
tracking technique measures an integrated change of frequency of the radio wave. This measuring process can be kept 
coherent for years as was demonstrated in missions Pioneer 10 and 11 
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The down-link of the Doppler tracking technique works as fohows. Spacecraft's receiver tracks the up-hnk carrier 
ft'equency using a phase lock loop. The signal received on-board of spacecraft with frequency 1^2 is coherently trans- 
mitted back to the ground DSN station and is received with a frequency v^. Calculating the down-link Doppler shift 
with the help of the light cone equation ([55]) . we derive 

^ 1 - fc32 ■ (/33 - fc32J^t3) , . 

V2 1 - fe32 • (/32 - fe32J^^2) ' ^ ' 

where the unit vector k-yi = r 32 1^32 is directed from the spacecraft to the Earth, and quantities with sub-index 3 refer 
to coordinates and velocity of the DSN station at the time of reception of the radio wave that was re-transmitted 
from the spacecraft toward the Earth at the time ^2- Technically, the procedure of re-transmission also includes 
multiplication of the re-transmitted frequency by a constant factor in order to prevent interference between incoming 
and outgoing radio signals [48|. We shall ignore this multiplicative factor for it changes frequency of the wave by a 
known number that is easily subtracted in the data processing but complicates theoretical description. 

Relationship between the emitted frequency vi and the observed frequency 1^3 is found from the condition that the 
phase of the radio wave in up- link is locked to that in the down-link. Hence, by applying (1^^ . we obtain 



^ yZ_V2_ ^ 1 - ^32 ■ (A3 - ^32^13) 1 - fc21 ■ (/32 " fc2lJ^^2) 
1^1 1 - fc32 • (^2 - k32^t2) 1 - fc21 • (/3l - k2l'Kti) 



(94) 



Because the light travel time between Earth and spacecraft is much smaller than the orbital period of Earth around 
Sun, we can neglect the orbital acceleration of the Earth and use the following equalities, fc32 — —^21 and (3^ — fii 
which are valid with a very good accuracy. In fact, the acceleration terms are taken into account in real data 
processing but they are not important in the present paper. We substitute the above equalities to ([94)) . and expand 
its denominator in Taylor series. It takes (|94l) to the following, simpler form 

— = 1 - 2fc2i • [132 - /3i - fc2l?^ (<3 -h)] , (95) 

where (32 — (3i is the relative radial velocity of spacecraft with respect to Earth. Introducing notation /32i = ^21 • 
{(32 — (3i) for the radial speed of spacecraft with respect to Earth, we write equation (|95|) as follows 

^ = l-2[/32i-J{(t3-ii)] . (96) 

Notice that (^21 > when spacecraft is receding from Earth, and (^21 < in case when it moves toward Earth. 

The up/down- link radio cycle is immediately repeated by transmitting to spacecraft at time of a radio wave 
with the same phase and frequency, 1^3, as the arrived radio wave. We notice that the technical procedure used in 
the DSN Doppler tracking is formally different as the signal is transmitted at with the reference frequency vi (not 
1^3) while the difference between i>3 and vi is time stamped and recorded [48]. Since the difference between i>3 and 
Vl is precisely known, we can always use in our theoretical development, frequency 1^3, as that transmitted to the 
spacecraft at time ^3. This simplifies theoretical understanding of the process of accumulation of the phase of a radio 
wave in the two-way Doppler tracking measurement and does not change the final formula for the integrated Doppler 
shift. 

The Doppler tracking process continues coherently following the same protocol. The wave, transmitted at time ^3, 
will be received at spacecraft at time with a frequency 1/4, re-transmitted back to Earth with the same frequency 
i/4, received at DSN station at time ^5 with frequency 1^5, transmitted back with frequency 1/5 to spacecraft, and so 
on. Thus, the wave with a reference frequency h'l emitted at time ti will be received at DSN station being shifted to 
frequency i'2n+i at time t2n+i after n up/down link bouncing cycles. Because we keep the phase of the wave locked, 
the relationship between the observed frequency, i^2ti+i, and the emitted frequency, vi, is given by the product of n 
factors 

V2n+1 _ ^2n+l l'2n-l V3 
V\ V2n-\ V2n-Z V\ 

where the corresponding calculation of the frequency ratios is done similar to that of v^jvx. Substituting these 
frequency ratios to (j97p . summing up and neglecting the acceleration terms, yield 

n 

= 1-25] \h^,2^-l - K {t2^+l ~~ t2.-l)] , (98) 
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where /32i,2i-i = ^21,21-1 ■ (/32i — /32i-i) is the relative radial velocity of spacecraft with respect to Earth at tinie t2i, 
fc2i,2i-i is the unit vector from the position r2i_i of the DSN station taken at time t2i-i along the light ray toward 
spacecraft's position at time t2i, /32i = 1^21/0, /32i-i = ^21-1/0, V2i is velocity of spacecraft at time t2i and V2i-i 
is velocity of the DSN station at time i2i-i- 

Let us recast equation ([98)l to the form that serves better to further discussion. We shall denote the observed 
frequency i^obs = i^2n+i, and equate the emitted frequency to the high-stabilized reference frequency vq of the 
master oscillator of atomic clocks, 1^0 = 1^1, driving the DSN radio transmitter. We shall also denote the overall time 
of the integrated Doppler measurement as t — t2i+i — ti. The Orbit Determination Program (ODP) adopted by 
NASA Jet Propulsion Laboratory (JPL) models (in the approximation under consideration) the integrated Doppler 
shift as follows 

J^mod = 1^0(^1- 2 J2 |32^,2^-?j ■ (99) 

It is more convenient for subsequent discussion to write (j99p in terms of an integral with respect to time by making 
use of the fact that the velocity of the DSN ground station and that of spacecraft do not change too much during 
the round trip of light from the DSN transmitter to the spacecraft and back. We can approximate the radial velocity 
with the increment 

^ = ^I^^2^At2. , (100) 
c 

where, amod — c{dl3/dt), is the modelled relative radial acceleration of spacecraft taken at time lying inside the 
time interval A.t2i — ^21-1 — ^22- Substituting (jlOOp to (|99|). and approximating the sum with an integral, we obtain 
exactly the equation used in ODP 
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t'mod = i^o 1 / a^od{t')dt' , (101) 



where t is the total span of observation, t' is the argument of integration and, for the sake of simplicity of notations, 
we had started the integration at the initial epoch, t = 0. We further notice that the accumulated phase of the 
radio wave (normalized to factor 2tt) emitted by the radio transmitter is 4>a — v^t, because the frequency of the 
master oscillator, is constant. On the other hand, the modeled total phase of the signal, is given by the integral, 
f^mod = /o Vmod{t')dt' , whcrc t'mod(i) is givcu by (|101l) . Performing integration by parts, yields 
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t 

{t^t')an,od{t')dt' . (102) 



Cosmological arguments tells us that the observed Doppler frequency must be described by equation (|98p which 
includes the term depending on the Hubble expansion as a linear function of time. This prediction differs from that 
obtained with the JPL's ODP formula ((99)) that was derived under assumption that the background spacetime is flat. 
The difference between the two equations can be written as follows, 

''"'^^ ~ ''"'"'^ = 2m , (103) 

where the term in the right side is obtained directly by summation of all contributions from time intervals in (|98p : 
Sr=i (^2i+i — ^2i-i) = 3^^- The phase difference between the observed and the ODP-modelled signals is obtained 
by integrating equation (|103p with respect to time 

Both equations, (jlOSp and (|104p . can be formally obtained from (|102p by replacing the ODP-modelled radial acceler- 
ation of spacecraft, Omod, with its cosmological counterpart, 

a{t) = a,„od(t) - cJf . (105) 

We notice that equation (|105p is a projection on the radius- vector connecting observer on the Earth with spacecraft. 
Therefore, the constant "anomalous acceleration" , Oanomaiy = — !Hc, of spacecraft should be interpreted as being 
directed toward observer but not to the Sun. 
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Equation (|103|) elucidates that the contribution of the Hubble expansion to the Doppler shift of frequency grows 
linearly with time t and, thus, can be measured after sufficiently long time interval if the phase of the radio wave 
used in the Doppler observation is kept locked during the time of observations (the coherent Doppler tracking). For 
example, during one year the fractional uncertainty between the observed and the ODP-modelled frequency as caused 
by the right side of equation (I103p . amounts to 2J{ x (1 yr) ~ 1.5 x 10^^". Long-term stability of atomic clocks 
is limited by a random walk frequency noise of Sv/i/ ~ 1.0 x 10~^^-\/t where time t is measured in days |5l!]. For 
the interval of time t = 1 year it is of ~ 2.0 x 10"^'' that is much better than the bias introduced by the Hubble 
expansion to the frequency drift of the coherent Doppler tracking. Analysis of the long-term stability of the coherent 
Doppler measurements provided by Anderson et al. [H, section VII F] yields a similar estimate for the error in 
measuring the left side of equation (|103p . After comparing the estimates for the left and right sides of equation (jlOSp 
we conclude that the Hubble drift of the radio frequency of the observed Doppler signal is detectable with the existing 
telecommunication technologies. 



D. Does the Pioneer effect measures the Hubble constant? 

When you have two competing theories that make 
exactly the same predictions, the simpler one is 
the better. 

Occam's Razor Prmctple 

Pioneer 10 and 11 were launched by NASA in the 1970s and radar-tracked for over 30 years. Experiment was 
carefully analyzed by J. Anderson with colleagues 0, |6l|. The data consistently indicated to the presence of an 
"anomalous acceleration" of spacecraft, that could be perfectly modeled by a "phenomenological" Doppler shift 
equation (see 0, Eq. 15]) 

= , (106) 

Vq c 

with ap = (8.74± 1.33) x 10~* cm/s^. It was noticed that the value of ap is very close to "Kc pointing out to possible 
explanation of the "anomalous acceleration" by the Hubble expansion. However, the sign of the effect (blue shift) 
was opposite to that having been naively expected from cosmology (red shift), and it was the main reason why all 
previous theoretical attempts to explain the blue shift of the Pioneer effect by cosmological expansion, have failed 
[lot [60j. Meanwhile, experimentalists were trying to explore more deeply the hypothesis that the anisotropic heat 
emitted by the spacecraft's generators is the primary cause of the effect [13, Ell i'^ spite of the original paper by 
Anderson et al. [2| had presented convincing arguments that the heat (and the other on-board generated systematics) 
cannot be responsible for the overall effect but only for its small fraction (~ 15% of the effect). Nevertheless, the 
efforts to explore the efficiency of the radiative force-recoil model continue. I am aware (private communication) of 
at least three works that are currently in progress, not as of this date announced or published, that in effect claim 
that the Pioneer anomaly can be explained in terms of heat radiation by the spacecraft. Definitely, the heat of the 
on-board generators should produce some effect on the orbital acceleration of the spacecraft. However, the present 
paper points out that there is a strong theoretical argument in favor that the prime cause of the Pioneer anomaly is 
deeply rooted in the fundamental gravitational physics of spacetime rather than in the thermal process. 

The theory of various time scales worked out in the present paper, is able to explain the Pioneer effect as caused 
by a secular drift of the frequency of a radio wave that is phase-locked in the DSN radio communication system 
observer-spacecraft. The drift of the radio wave frequency is caused by the cosmological expansion that affects the 
rate of electromagnetic time scale. A, as compared with that of the atomic and ephemeris time scales, t, (see equation 
dnil)). It explains both the magnitude of the frequency drift and the blue-shift sign of its time derivative. The effect 
is absent in Minkowski spacetime where all time scales (atomic, ephemeris, and electromagnetic) measured by a set 
of fundamental observers, have the same rate. 

It is worthwhile to notice that 0] attempted to fit the Pioneer effect by introducing a quadratic time divergence 
of atomic clocks from the ephemeris time of gravitational equations of motion. This did not work out and this 
negative result is in a full accordance with our prediction that the rate of atomic and ephemeris time scales have no 
divergence in expanding universe. On the other hand, Anderson et al. [2j were successful to fit the Pioneer effect 
with the phenomenological quadratic-in-time model that was applied only in equations of light propagation. However, 
neither at that time nor later the quadratic-in-time model did not find theoretical justification and remained purely 
phenomenological toy-model [37| . Our theory provides a solid theoretical framework supporting the quadratic-in- 
time model originating from the non- uniformity of the electromagnetic time scale, \ = t -\- \/2'Kt^, realized by a 
photon bouncing between two 'radio mirrors' (observer-spacecraft) in an expanding universe. This non-uniformity of 
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electromagnetic time scale leads unequivocally to the Doppler shift equation (I103P and provides theoretical explanation 
to the experimentally-established equation (I106|) . 

Equation (I103P can be also reformulated in the form of (jlOSp . and interpreted in terms of the constant "anomalous 
acceleration" , aanomaiy = — 3^c, of spacecraft directed toward observer. If spacecraft is moving away from observer 
the "anomalous acceleration" leads to a tiny blue shift of the Doppler frequency of the spacecraft's signal on the top 
of a large red shift caused by the outgoing motion of spacecraft. We also notice that equation (|105|) is a projection 
on the line connecting observer with spacecraft. Therefore, the "anomalous acceleration" is directed toward observer 
on the Earth but not to the Sun. It agrees with the observed direction of the "anomalous acceleration" . Equation 
(|105|) also tells us that the "anomalous acceleration" is constant everywhere within the solar system and far outside 
of its boundaries. This is because the "anomalous acceleration" is not a physical force but the effect of the cosmic 
expansion which is homogeneous and isotropic everywhere in the FLRW universe. 

The "anomalous acceleration" is the most common phenomenological interpretation of the Pioneer effect adopted 
in most literature [131 ■ emphasize that our theory associates the physical origin of the Pioneer effect neither 
with the "anomalous acceleration" nor with an "anomalous gravitational force" but with the different rates of the 
time scales caused by the cosmological expansion. The search for the "anomalous gravitational force" was a matter 
of many theoretical speculations aimed to explain the Pioneer effect with a "new physics beyond Einstein" [37| . 
Our standpoint is that the Pioneer effect can be explained within the framework of general relativity as a natural 
consequence of the (locally-observed) Hubble expansion without resorting to alternative theory of gravities or other 
exotic phenomena. 

Comparing (|106|) with (|103p yields a theoretical prediction for ap = "Kc. Taking the WMAP value for "K [2^ leads to 
a^MAr ^ (6.9 ± 0.3) X 10"** cm/s^ On the other hand, substituting the numerical value for ap = (8.74 ± 1.33) x IQ-^ 
cm/s^, measured in Pioneer experiment 0, yields a corresponding value of the Hubble constant M = 90.0 ± 13.7 
km/s/Mpc that is larger than the value of J£ = 71.0 km/s/Mpc derived from WMAP data. More recently obtained 
value of ap = (8.17 ± 0.02) x 10"* sm/s^ ^ yields, J{ ^ 84.2 ± 0.2 km/s/Mpc, that still noticeably exceeds the 
WMAP value. The difference is about 20% and may be explained, indeed, by accounting for the residual thermal 
emission of spacecraft's generator. 

Cassini spacecraft was also equipped with a coherent Doppler tracking system and it might be tempting to use 
the Cassini telemetry to measure the "anomalous acceleration effect" which is predicted by the same theoretical 
equation (|105l) and, if WMAP value for 3-C is implemented, has the magnitude of ~ 7 cm/s^. Unfortunately, there are 
large thermal and outgassing effects on Cassini that would make it difficult or impossible to say anything about the 
"Pioneer anomaly" from Cassini data, during its cruise phase between Earth and Saturn Q. Due to the presence of 
the Cassini on-board generated systematics, the recent study (23{ of radioscience simulations in general relativity and 
in alternative theories of gravity is consistent with a non-detection of the "anomalous acceleration" effect. 



X. THE CONCLUDING REMARKS 



A. Cosmological Effects in Local Inertial Coordinates 

There are more things in heaven and earth, 
Horatio, than are dreamt of in your philosophy. 

Shakespeare. Act I, scene v, Hamlet 

The reader may be wondering how it is possible that the Hubble expansion can be measured locally. It seems 
to be in contradiction to a common wisdom which states that any curved manifold admits a construction of local 
inertial coordinates along any time-like worldline such that all effects produced by the manifold's curvature vanish at 
the origin of the coordinate chart. This is because the Christoffel symbols are nil at the origin of the local inertial 
coordinates [47j . 

The present paper does not argue against this well-established mathematical fact. The local inertial coordinates 
in cosmolo gica l spacetime can be built in accordance with the mathematical principles of the Riemannian geometry 
pH [2^ [25I |30( . The question is not about mathematics but about how the coordinate description of motion of test 
particles and light relates to physical observables - the proper time and frequency of light - measured by an observer at 
the origin of the local inertial coordinates. Physical treatment of this problem in the linearized Hubble approximation 
is given below. Many technical complications will be passed over. More specifically, we shall drop out terms of the 
order of Jf^ from equations. We shall also neglect the effects of the gravitational field of the localized astronomical 
system which are incorporated to the metric tensor perturbations laj3 ■ They can be easily added up by making use of 
the principle of linear superposition as confirmed by calculations given in sections fVl IVII of the present paper. 
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We start from the FLRW metric (|6]) and make transformation ((59)) of spatial coordinates, ~ R{t)x'^. We obtain 
the FLRW metric in the new coordinates, 

ds^ = -c^df - 23<r'dtdr' + Sijdr'dr^ , (107) 

where all terms of the order of Jf^ have been omitted. We remind that the cosmological time t is the proper time of 
the Hubble observers having fixed spatial coordinates 2;* = const.. Clocks of the Hubble observers are syiichronized 
with the help of Einstein synchronization procedure based on exchange of electromagnetic signals [H, [g^] , and the 
rate of the clocks is uniform. 

Strictly speaking, coordinates {t, r*} are not locally inertial because not all ChristofFel symbols in these coordinates 
vanish. In particular, the only non- vanishing ChristofFel symbol is T'^ij = {Ji/cjSij. This Christoffel symbol causes a 
non-incrtial force and is responsible for the accelerated motion of photons in coordinates {t, r*} as shown in equation 
of motion ([5^ . On the other hand, the non- vanishing ChristofFel symbol affects equations of motion of slowly moving 
test particles only in the post-Newtonian terms of the order of Jf/c^ which are negligibly small in the solar system. 
This is the reason why the coordinates {t, r*} can be treated as inertial in the Newtonian limit of the slowly moving 
particles. It corresponds to the discussion given in section IVII Al 

Metric p07p can be further simplified by making use of time transformation, 

A = < + if o,,rV^ + . . . , (108) 



where A is a new coordinate time, and the dots denote residual terms of the higher order in the spatial coordinates 



and the Hubble parameter. After substituting the above-given value for = CK/cjSij, the explicit form of the 
coordinate transformation becomes 



A = <+^r2, (109) 

where the residual terms have been dropped out. The reader should notice that equation ()109|) apparently indicates 
that the spatial hypersurfaces of constant time, A, do not coincide with the hypersurfaces of simultaneity of the proper 
time t of the Hubble observers. 

The metric (|107p has in the new coordinates {A,r*}, apparently Minkowskian form 

ds^ = -c^dX^ + S.jdr'dr^ . (110) 

It reveals that coordinates {A,r*} are locally inertial coordinates introduced on the cosmological FLRW spacetime 
manifold in the neighborhood of the worldline of a Hubble observer having fixed spatial coordinates = 0. All 
ChristofFel symbols disappear at the origin of the local inertial coordinates {A,r*} that can make an impression that 
all cosmological effects which might be measured are shifted to the residual terms of the order of Jf ^ . This conclusion 
is invalid. 

The fact is that the coordinate time A is not the proper time t of the Hubble observer outside of the observer's 
worldline, and is not directly measurable. Therefore, the coordinate description of motion of test particles with respect 
to the coordinate chart, {A,?"*}, is not given in terms of physically-measurable time, t, and should not be interpreted 
like the motion of particles in Minkowskian spacetime. To make this important point more clear, let us consider a test 
particle moving freely from the origin of the local inertial coordinates. Equation of motion of the particle is geodesic, 
which is given in the metric (IllOp by the first Newton's law 

= 0. (Ill) 



dA2 

This equation can be integrated easily 

r' = v'X , (112) 

where = dr"^ /d\ is a constant velocity of the particle. We have to replace the unobservable coordinate time. A, with 
the directly-measurable proper time t of the observer which is located at the origin of the coordinates. Relationship 
between A and t is obtained from the coordinate transformation (I109P which must be performed at the current position 

A = t + ij{^. (113) 



of the particle given by equation (jll2 
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Solving by iterations, we get 




(114) 



where the dots denote residual terms of the higher order in t. Equation (|114p demonstrates that the coordinate 
time, A, diverges quadratically from the proper time, of the observer as the particle moves away from the origin 
of the local coordinates. The divergence between the two time scales, A and t, is a reflection of the symmetry of the 
underlying cosmological manifold pointing out that the coordinate simultaneity in the local inertial coordinates does 
not correspond to the physical simultaneity of clocks of the Hubble observers in the expanding universe. 

The divergence between A and t is negligibly-small for slowly-moving particles due to the presence of the post- 
Newtonian parameter, (w/c)^ ^ 1, in combination with the Hubble parameter IK. However, if the particle is a photon 
(radio wave), its velocity v = c, and equation (jll4l) becomes 



This equation is identical with (1551) and elucidates our derivation of the propagation of light given in section IVHBI 
We conclude that the Hubble parameter enters equations of motion of test particles in the local inertial coordinates in 
cosmology. Time coordinate, A, is a convenient mathematical parameter along the path of test particles but it must 
be expressed in terms of the observable proper time t of the Hubble observer in order to make physical sense of the 
parametrization. 



Pioneer effect caused a keen interest in scientific community of experimental and gravitational physicists. There were 
several proposals submitted to NASA and ESA [12, 61] describing dedicated space missions to test the "anomalous 
force" and "celestial anomalies" in the solar system. All proposed missions requested the launching spacecraft to deep 
space beyond the orbit of Jupiter. 

According to our theory, the testing of the "anomalous force" has nothing to do with the violation of general 
relativity. The effect of the "anomalous force" detected in Pioneer 10 and 11 missions, is associated with the non- 
uniformity of electromagnetic time scale (Marzke- Wheeler clock [i^) as compared with atomic time scale (atomic 
clock) and the ephemeris time scale (gravitational clock). Therefore, the effect can be tested in any space mission 
having highly stabilized atomic clock which long-term stability is measured with a coherent Doppler tracking system. 
The most restrictive requirements to the mission are: 

1. the suppression and/or precise calibration of the on-board systematics like heat, gas leakage, etc., 

2. the reduction of the impact (or precise measurement with on-board detectors) of the non-gravitational forces 
caused by the solar wind, light pressure, dust, etc., 

3. deployment of ranging measurement technique that ensures an independent determination of the space probe 
position and velocity. 

PHARAO project may be an excellent candidate for measuring the effect of the local cosmological expansion which 
has a well-predicted magnitude, sign and direction in accordance with our equation (I103|) . Another possibility is to 
deploy a radio transmitter in a package with laser retro-reflector on the visible face of the Moon in order to conduct 
the coherent Doppler and laser ranging measurements to test (llOSp . It can be also achieved with the high-precision 
Doppler and satellite laser ranging of a telecommunication spacecraft launched at a geostationary orbit. Constellation 
of GPS satellites may be useful as well but this should be explored more deeply as the GPS satellite orbit degrades 
as time goes on. 

It may happen that measuring the local value of the Hubble constant can be achieved with the optical flywheel 
oscillator [4g. The oscillator represents an optical resonator whose current instability is ~ 3.5 x 10~^^t~^/^ . The 
frequency of this type of oscillator, as compared with that of atomic clocks, is expected to drift due to the Hubble 
expansion as J-Ct in accordance with (jl03p . In principle it can be measured already after 1000 seconds. However, 
the theory of the effect of the Hubble local expansion in optical oscillators should be worked out in addition to the 
calculations of the present paper. Moreover, frequency instability of optical oscillator is caused by many other factors 
which should be carefully studied and subtracted. 
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